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The holonomy correction is one of the main terms arising when implementing loop quantum 
gravity ideas at an effective level in cosmology. The recent construction of an anomaly free algebra 
has shown that the formalism used, up to now, to derive the primordial spectrum of fluctuations 
was not correct. This article aims at computing the tensor spectrum in a fully consistent way within 
this deformed and closed algebra. 
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II. THEORETICAL FRAMEWORK 
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Nonperturbatively quantizing General Relativity (GR) 
in a background-invariant way is obviously an outstand- 
ing open problem of theoretical physics. Loop Quantum 
Gravity (LQG)Js a promising framework to perform this 
program (see [l| for introductory reviews). It has now 
been demonstrated that different approaches, based ei- 
ther on quantizations (covariant or canonical) of GR, or 
on a formal quantization of geometry lead to the very 
same LQG theory. Experimental tests are, however, still 
missing. Trying to find possible observational signatures 
is a key challenge and cosmological footprints are known 
for being one of the only possible paths toward a real 
experimental test of LQG. It is very hard to make clear 
predictions in Loop Quantum Cosmology (LQC) using 
the full "mother" LQG theory. General introductions to 
LQC can be found in 0. This study focuses on an effec- 
tive treatment taking into account recent results on the 
correct algebra of constraints. We first review the theo- 
retical framework. The spectrum is then derived. Some 
conclusions and consequences are finally underlined. 
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One of the fundamental quantum corrections expected 
from the Hamiltonian of LQG arises from the fact that 
loop quantization is based on holonomies, i.e. exponen- 
tials of the connection, rather than direct connection 
components. Based on a canonical approach, the theory 
uses Ashtekar variables, namely SU(2) valued connec- 
tions and conjugate densitized triads. The quantization 
is obtained through holonomies of the connections and 
fluxes of the densitized triads. This is the key ingredient 
of the effective approach. The cosmological equations are 
modified so as to account for the loop basis of the theory. 

The main consequence of the holonomy correction 
on the cosmological background is to induce a bounce. 
The evolution is not singular anymore and the Big Bang 
is replaced by a Big Bounce. The next step consists 
in studying the propagation of perturbations within 
this modified background. In cosmology, perturbations 
are of three different types : scalar, vector and tensor. 
We focus here on the tensor modes that are directly 
gauge- invariant. Quite a lot of works have already been 
devoted to tensor modes in this framework [3|. Beyond, 
the phenomenology of LQG is now a well established 
field (see [|[ for a review). Unfortunately, a recent study 
[EJ has shown that the previously derived spectra are 
most probably incorrect. 

The key issue relies in the closure of the algebra of 
constraints. Due to general covariance, the canonical 
Hamiltonian is a combination of constraints C/. Con- 
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sistency requires that the constraints are preserved un- 
der the evolution they generate. This is ensured in the 
classical theory by the closure of the Poisson algebra of 
constraints 



{C I ,Cj} = f K IJ (Al,E9)C K , 



(1) 



where Ci, 7 = 1,2, 3, are the Gauss, diffeomorphism and 
Hamiltonian constraints and f K IJ (A l b ,Ef) are structure 
functions which, in general, depend on the phase space 
(Ashtekar) variables (A 3 b ,Ef). They form a first class 
set. Otherwise stated, the gauge transformations and 
evolution generated by the constraints define vector fields 
which are tangent to the sub-manifold defined by the 
vanishing of constraints. 

In LQC, quantum corrections are introduced as effec- 
tive modifications of the Hamiltonian constraint. This 
generates anomalies: the modified constraints C/ do not 
form a closed algebra anymore: 



{C? ) C*} = f K IJ (Ai,E?)C% + A 



i j- 



(2) 



The anomalous terms Au are removed by carefully ad- 
justing the form of the quantum correction to the Hamil- 
tonian constraint through the addition of suitable "coun- 
terterms" that vanish in the classical limit. This has been 
done in Q , following the approach of [f| . 

In the classical case, the Poisson brackets between the 
constraints read as: 

{D M [N?],D {mii)) [N$}} = 0, (3) 
{H {mh)) [N] , D {mM)) [N a ]} = -H {m ^ [6N a d a 6N] , (4) 



{H(m+g) Wl] , -^(m+g) [-^2] } = D { , 



in+cj) 



fd a (SN 2 



6N1) , (5) 



where (rn+g) stands for gravity and matter. The quantum 
corrections are included at the effective level by replacing, 
as usual, in the Hamiltonian constraint 



sm(p'yk) 
PI 



(6) 



The important result of Q is that the quantum-corrected 
algebra is described by a single modification: 



{H(m+g)Wl], Hf jn ^ j [N 2 ]} = J~2 D( m+g ) 
where 
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d a (5N 2 - 5N X ) 



fl = cos(2/l7fc) = 1-2— 

Pc 



(7) 
(8) 



The fl factor encodes the quantum correction, k being 
the homogeneous Ashtekar connection and p, being pro- 
portional to the ratio between the Planck length and the 
scale factor. The Mukhanov-Sasaki Q equation of mo- 
tion for gauge-invariant perturbations of scalar and ten- 
sor types vq(t) can be explicitly derived. In conformal 
time, the propagation of tensor modes is given by 



u T -nw 2 v T 



-vt = ; zt 



a 
fl 



(9) 
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FIG. 1: Evolution of fl and its derivative with respect to 
conformal time. The density where fl vanishes is half the 
critic density whereas fl' vanishes at the bounce. 



where prime means differentiation with respect to confor- 
mal time. This leads to the following equation of motion 
for tensor perturbations, defined via vt = zt x h l a : 



!>:: - /v, ( ■!.// - ~ I -nvx = 0. (10) 

where :— a' /a is the conformal Hubble parameter. 



III. POWER SPECTRUM 

This equation being known, it is possible to investigate 
the associated primordial power spectrum. This is the 
fundamental ingredient for phenomenology. The back- 
ground dynamics is not modified by the fl term. How- 
ever, the perturbations will of course undergo a different 
evolution. 

We use the Fourier transformed version of Eq. (fTU|) : 



h' + flk 2 h = 0, 



(11) 



where the indices have been skipped for simplicity. The 
behavior of fl and fl' is displayed in Fig. [TJ One can 
immediately see that fl vanishes for p = Pc/2, where 



Pc = 



V3 



32tt 2 7 



0.82m 



pi- 



(12) 



In addition, fl becomes negative- valued, leading to an 
effective change of signature of the metric (Euclidean 
phase) around the bounce. Equation (|11[) is apparently 
ill-defined as fl' /fl — > oo at r) = rf' and -q = rf + \ the 
values of conformal time when p = p c j1 before and after 
the bounce, respectively. However, regular solutions do 
exist by rewriting Eq. (|lll) as: 



ti = fig ; g' = -2,J4°g - k 2 h, 



(13) 



which is regular. 
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FIG. 2: Evolution of the scalar field (upper panel) and the 
scale factor (lower panel) as a function of cosmic time (the 
bounce corresponds to t = 22693). The parameters are m = 
10" 3 Mpi and x B = -1.5 x 1(T 3 . 



The same set of equations in cosmic time are: 

h= —g ; g = -2Hg h, 

a a 



(14) 



where dot means differentiation with respect to cosmic 
time and H is the usual Hubble parameter. The dynam- 
ics can also be recast in a single second order equation: 
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23fg' + (2JT' + flk 2 )g = 0. 



(15) 



Whatever the chosen form, either (flUf . (fT4")l or (TT5"1) . 
the evolution can be computed numerically. Of course, 
the propagation of modes has to be coupled with the 
background evolution which is drastically modified by 
the holonomy corrections that are at the origin of the 
bounce. The cosmological background evolution is ba- 
sically driven by a single scalar massive matter field of 
mass to. We define 



x :- 



m(f> 



and y := 



(16) 



which respectively represent the density of potential and 
kinetic energy normalized so that x 2 B + y B = 1 at the 
bounce. The free parameters of the study are therefore 



7f 



|hf 
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FIG. 3: Mode amplitudes as a function of time, corre- 
sponding (from top to bottom at i = 20000) to k = 
10 2 , 10 1 ' 5 , 10\ 10 ' 5 , 1, 10 -0 ' 5 , 10 _1 , 10" 1 ' 5 , 10" 2 , 10" 2 ' 5 , 10 -3 . 
The parameters are m = 10 _3 Mpi and xb = —1.5 x 10 -3 . It 
should be noticed that the initial conditions fore each mode 
are specified long before the time interval of this plot. 



to, xb (the value of x at the bouce) and the relative sign 
of xb and i/b- Interestingly, if the initial conditions for 
the background are specified at any time, long enough 
before the bounce, the probability of \xb\ is strongly 
peeked around a given value of order m (in Planck units), 
with sign(a;s) = sign(ys) (the detailed probability dis- 
tribution for xb will be studied somewhere else [8]). For 
numerical reasons it is better to specify computational 
initial conditions for the background before the bounce 
rather than at the bounce. Because of the peaked proba- 
bility, the resulting xb is always close to the same value. 

It is also necessary to assign a numerical value to 
the scale factor, a at some point. This choice has of 
course no physical consequences but has to be taken 
into account for the interpretation of the meaning of 
the wave vectors k 7 since they are expressed in the 
coordinate space and not in the physical space. The 
explicit choice made was a = 1 at the bounce, which is 
numerically easier than the usual normalization at the 
nowadays value. 

In Fig. [21 the evolution of the scalar field and scale fac- 
tor are shown for some typical parameters. As expected, 
the oscillations of the scalar field are amplified before the 
bounce, because the negative Hubble parameter acts as 
an anti-friction term. Then, just after the bounce, the 
Hubble parameters becomes positive and large, acting as 
a huge friction and therefore leading to slow roll inflation. 

The amplitudes of some Fourier modes of h are plotted 
in Fig. [3] They are obtained by chosing the Minkowski 
vacuum as the initial state, since z" / z — ¥ in the remote 
past. 
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Before the bounce, for k 2 > z"/z, \h\ 2 = \/{2ka 2 ). 
When z" I z s» k 2 or z"/z > k 2 , \h\ 2 grows quicker. 
Since the amplitudes of smaller k start growing quicker 
before the amplitudes of larger k, this adds up to a 
collecting effect that brings all modes up to a certain 
k m max t<fB (\/z"/z) up to the same amplitude. Af- 
ter the bounce, the amplitudes oscillate until k 2 ^> z" / z 
when we get v oc a (as can bee seen from Eq. (|9])) and 
therefore h = constant. 

Finally, the power spectra for different cases are pre- 
sented in Fig. 2J The main features are the following: 

• a flat (scale invariant) infrared limit, 

• an oscillating intermediary part, 

• an exponential behavior in the ultraviolet limit 
(starting around k = 2 independantly of m) . 

This obviously exhibits important deviations, with re- 
spect both to the standard GR case and with respect to 
previous LQC computations without the term. Al- 
though surprising at first sight, the exponential diver- 
gence in the UV limit might not be catastrophic as 
physics at very small scale is anyway not described by 
the primordial power spectrum. 

Furthermore, this ultraviolet behavior can be checked 
analytically. In the large k limit of Eq. (j^J), the 
WKB conditions are satisfied in the euclidean phase 
around the bounce. More precisely, those WKB con- 
ditions are met for r\ € \rf~^ + el \t)^ — e fc + ^] with 

4 ±) ~ ( fc2 = ?7 (±) |)~ 1/3 - The Mukhanov-Sasaki 
function can be approximated by 

v T = v+ e tk f^ d ^ + v-e^ k ^ VUdv . (17) 

As O is negative-valued during the euclidean phase, the 
tensor mode is dominated by its exponentially-growing 
solution 

( r (+) -^ +) ,— \ 

h oc exp [k / VWldrj . (18) 



This can also be seen in Fig. [3] where the amplitude of 
large k modes grows rapidly in the vicinity of the bounce, 
where O < 0. 



IV. DICUSSION 

This study implements in a consistent way the mod- 
ified algebra induced by holonomy corrections in the 
calculation of the primordial tensor power spectrum. 
Thanks to numerical calculations, it was possible to solve 
the equation of motion for gravitational waves. The re- 
sulting spectrum exhibits specific features. Of course, 
this raises important questions. First, the well known 
problem of trans-planckian modes in inflation (see, e.g., 
9]) should be treated with a specific care in LQG where 
the very meaning of a length smaller than the Planck 
length is dubious. If the number of c-folds of inflation is 
chosen (by appropriately setting a very small fraction of 
potential energy density at the bounce) to be just above 
the minimum required value, then modes relevant for 
the CMB are still sub-planckian and the approach makes 
sense anyway. In other cases, the effective theory might 
just breakdown. With the normalization chosen in this 
work the trans-planckian window corresponds to k > 1. 
Second, the propagation of modes through the euclidean 
phase is not straitghforward [Io| • Strictly speaking, there 
is no "time" in that region and the concept of evolution is 
not well defined. In this work, we have deliberately cho- 
sen to withdraw the conceptual issues associated with 
the transition between hyperbolic and elliptic solutions 
and to focus on a well defined mathematical solution. An 
alternative approach, based on the BKL conjecture, will 
be studied later [ll|. An analogous study should also 
be performed for scalar modes. The regularization trick 
used here, however, does not apply directly, and other 
methods have to be constructed. We stress out that the 
case of scalar modes with holonomy corrections has been 
studied in [l2| and 

El 

but in different settings for the 
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background ; for the study of 1121 is restricted to super- NSF grant PHY- 1205388. 
inflation while the study of [13[ considered a dust-like 
bouncing Universe. Finally, those results will have to be 
compared with forthcoming studies based on other very 
recent approaches to LQC [l4{. 
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